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We propose a new selection criteria for predicting the most probable wavefunction of the universe 
that propagates on the string landscape background, by studying its dynamics from a quantum 
cosmology view. Previously we applied this proposal to the SUSY sector of the landscape. In this 
work the dynamic selection criterion is applied to the investigation of the non-SUSY sector. In the 
absence of detailed information about its structure, it is assumed that this sector has a stochastic 
distribution of vacua energies. The calculation of a distribution probability for the cosmological 
constants A e //, obtained from the density of states p, indicates that the most probable wavefunction 
is peaked around universes with zero A e / / • In contrast to the extended wavefunction solutions found 
for the SUSY sector with N- vacua and peaked around A e // — -^r, wavefunctions residing on the 
non-SUSY sector exhibit Anderson localization.Although minisuperspace is a limited approach it 
presently provides a dynamical quantum selection rule for the most probable vacua solution from 
the landscape. 

PACS numbers: 98.80.Qc, 11.25.Wx 



Recently we proposed to use quantum cosmology 
on the background of the landscape of string theory 
in order to find the most probable wavefunction of 
the uni verse 0. This proposal provides a new selec- 
tion criterion for predicting the most probable wave- 
function of the universe based on the dynamics of the 
superspace. In |l| we focused our attention on the 
minisuperspace of the supersymmetric (SUSY) vacua 
sector of the landscape. The SUSY and non-SUSY 
sectors of the landscape were assumed disconnected. 

Here we investigate solutions for the most probable 
wavefunction of the universe by extending our previ- 
ous analysis to the non-SUSY sector of the landscape. 
The minisuperspace approach is a reasonable approxi- 
mation for both SUSY and non-SUSY sectors because, 
as shown below, the previous assumption that the two 
sectors are disconnected is fully justified. The purpose 
of this work is to examine the question: Can we make 
a prediction for the most probable and stable solution 
for the wavefunction of the universe propagating on 
the background of the non-SUSY sector of the land- 
scape? 

Preliminaries and Setup: Consider the array of 
landscape vacua solutions in the non-SUSY sector to 
be a crystal. The array of vacua are parametrized by 
the potential V{<f) of the moduli field which collec- 
tively characterizes all the relevant degrees of freedom 
of string fluxes and fields 4> = [^n]- The vacua are 
located at fa, i = 1,...,N. The energies of the non- 
SUSY vacua, \((f>i) = Xi, can take any value in the 
range Aj € [0, ±W] where W — 0{M^) is an energy 
scale related to the Planck or string constant scale. 
The crucial assumption is that the vacua energies on 
the non-SUSY sector are stochastic, namely Aj's are 
drawn randomly in the interval [0, ±W]. The spacing 
between different vacua could also be a random 'sprin- 
kling' of sites. In such a setup, the non-SUSY sector 
of vacua is a randomly disordered lattice of configu- 
ration (super-)space of moduli <j), which we name the 



'superlattice'. The goal is to use quantum cosmology 
on the landscape background for dynamically select- 
ing the most probable wavefunction of the universe 
propagating on the 'disordered superlattice'. 

Wavefunction of the Universe on the Stochastic 
Minisuperspace of the non-SU SY Sector: Following 
the proposal of0 let us define the non-SU SY minisu- 
perspace to be restricted to the 'disordered superlat- 
tice' of the stochastic non-SU SY sector V(4>) of vacua 
4> and to homogenous flat 3— geometries 



ds 2 = [-Ndt 2 + a 2 {t)dx 2 } , 



(1) 



with M being the lapse function. The non-SU SY min- 
isuperspace is parametrized by the variables, [a, <f>]. 
The Lagrangian for the system with variables [a, <p] 
receiving contributions from both: gravity (L g ) and, 
moduli 'superlattice' (L^), is L = L g + where 
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\I/[a, <j>] denotes the wavefunctional of the universe 
propagating on the minisuperspace background [a, </>]. 
The Hamiltonian constraint on the wavefunctional is 
obtained by varying the combined action with re- 
spect to the lapse function M.\n the usual manner p^jj 
HE E3- promoting p a , pa, into operators, gives the 
Wheeler-De Witt (WDW)ja,H2 equation 
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where a is replaced by a = e a . (A treatment of the 
subtleties related to the normal ordering of opera- 
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tors and cross-terms in the Wheeler-DeWitt equation, 
Eqn. (3), can be found in Ref. [28,29,30]) 

Rescaling <j> to x = e 3a cf>, (and the other relevant 
quantities in the potential), in order to separate vari- 
ables in @, $f[a, <j>] can be decomposed in modes 



¥(a,x) = '£ J 1> j (x)F J (a). 
Replacing @ into © and using 

7i{x)%jjj{x) = ijtpj(x) where 
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results in 
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where 'hat' denotes the rescaled ij — e ea €j and from 

hereon-^ 2 - is absorbed into ej given in fundamental 
units. The a equation of motion is obtained by vary- 
ing the action, S resulting from Eqn. @ with respect 
to a 

a+l [a 2 + ( x 2 -V(x))e- 6a ] = 0. (7) 

It is easy to check for consistency that the a equation 
of motion is the Friedman equation for the expansion 
The 'eigenvalues' £j are obtained by solving the 
'Schrodinger' type equation, Eqn.( [5J, for the held 
tpj(x) propagating on the disordered superlattice V(x) 
with N lattice sites Xi and stochastic vacua ener- 
gies Xi. The hamiltonian H{x) = Hq + Hi con- 

tains two pieces: the diagonal part H — — Vq 
where Vo(xi) = A; and the short range interac- 
tion between neighbor vacua [xi,Xj], (the nondiag- 
onal terms), Hj = Vj .Therefore the potential V(x) 
includes two terms V(x) = Vq + Vi(x). The short 
range interaction allows spreading of the wavefunc- 
tion. Let us assume the nearest neighbor approxi- 
mation for the short range interaction, for example 
think of tuneling to the nearest neighbors for simplic- 
ity, and introduce a (dimensionless) (5-correlated white 
noise for the interaction term, i.e < Vi >= and 
< Vi(xi)Vi(xj) >= TS(xi — Xj) where < ... > denotes 
ensemble averaging with respect to [xj]. (Note that 
we do not assume < Vq >= A to be centered around 
zero). Due to stochasticity the potential and eigen- 
values can not be written in an exact form. In most 
cases, the details of V(x) = Vq + Vr do not matter 
since all the results are obtained only in a probabilis- 
tic manner. This is a complicated N— body problem 
because of the multiscattering of tpj( x ) among many 
sites of the 'superlattice' but relevant quantities are 
calculated from probability distributions which can be 
found exactly. 



There is a vast amount of literature on the prop- 
agation of wavefunctions on disordered lattices and 
the treatment of solutions for types of Eqn.© can 
be found in many papers, 0, 0, IE IE 0- P- Ander- 
son was awarded the Nobel prize for his investiga- 
tion into this issue. Independent of the specifics of 
the model, the conclusion is that the wavefunctions 
with energies below a certain scale set by the disorder 
strength do not propagate on random lattices with 
short range interactions. Instead, the wavefunction 
soon gets localized around a lattice site Xj with vac- 
uum energy Xj, a phenomenon known as Anderson 
localization[3| .Localization is purely a quantum me- 
chanical effect and it occurs because of the destruc- 
tive interference of the phases of the wavefunction 
from multiple scattering among sites. As shown in 
[3], this phenomenon always occurs when: disorder 
is present and, the interaction term Vi is sufficiently 
short-range, for all levels with energy below the dis- 
order strength set by F. The ID and 2D disordered 
systems are critical in the sense that localization oc- 
curs for all energy levels, independent of the disorder 
strength Y. Many examples of this phenomenon have 
been succesfully applied to lattice QCD0, observed 
experimentally in condensed matter systems and 
later derived from the (RMT) Wigner-Dyson theory 
of random matrices 0, H II S llS E3 > with the NXN 
matrices obtained over many realizations of the ran- 
dom potential^!]]. Unfortunately there is a no one-to- 
one correspondence between the localization site and 
the energy, since vacuum energies are stochastic vari- 
ables with values assigned randomly from the interval 
[0, ±W], and W = O(Mp) for our case. However, con- 
clusions can be drawn about the energy of the local- 
ized wavefunctions from their distribution probability 
shown below. 

Let us consider the distribution function of the 
vacua energies from the interval [0, ±W] by P(X). 
P(X) and therefore P(H(x)) can be a Gaussian dis- 

tribution P(X) = 7= e «r with ensemble aver- 

aged mean values < l|Aj|AT >= A and width Y =< 
l\Yi\N >, or; when disorder is large Y = 0(W) by 
a flat distribution P(X) = i. Corrections to the 
unperturbed energy for the wavefunction localized 
around Xj, along with the evaluated Green's function 
< j\G\j >=< — Hx)^ 1 ^ > can be estimated 
by the usual perturbation theory |3J for weak disor- 
der and by Wigner-Dyson RMT methods for the case 
when disorder is large and perturbation theory breaks 
down. The latter is a more elegant and transparent 
method of calculation since averaging < ... > is done 
by integrating with respect to H (x) with probability 

weig ht P[H{x)] ~ e— rather than ensemble 
averaging over vacua [xj]. The case of lattices with 
stochastic energy distributions can be found in many 
reviews 0, 0, IE B thus below we outline only the 
main steps of the calculation relevant to our purposes. 
Let us denote the shift of the unperturbed energy at 
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'site' Xj by Sj = <5j + ij. Perturbed energy 6j and 
Green's function are 

ej w Aj - Sj - ij 
G-ljj = (e-e j ) 

(8) 

with Sj the usual second-order correction that for sim- 
plicity we may take to vanish. With perturbation the- 
ory methods i'jj — 2-^ = iHk^j—j 2 —, where Vj jh 
denotes interaction between sites Xj and Xk- The 
same results are obtained by calculating Sj with the 
Wigner-Dyson RMT formalism, (for techincal reasons 
up to self-energy corrections in most cases), 7j ~ 2ze£. 
Assuming that, independently of the 'site' Xj (white 
noise), the second correlation probability momentum 
falls rapidly fast away from Xj then the most prob- 
able and averaged mean value for lj are roughly 
the same. Thus we replaced it with its mean value 
T =< 1 1 r j- 1 iV >, in the expression < V(xj)V(x k ) >« 
T5(xj - x k ).). 

The eigenfunction tpj (x) localized around Xj is given 

by 

H 

Due to Anderson localization if)j(x) can not propa- 
gate from the non-SUSY sector to other sectors of 
the landscape. (Therefore the assumption that the 
two sectors are disconnected seems fully justified). Lo- 
calization is a demonstration of purely quantum me- 
chanical effects. Classically the amplitudes of scatter- 
ing/tunneling through many vacua would be multi- 
plied over all the 'sites' 1..N which means that the 
final transition amplitude for propagation through a 
periodic lattice would be zero. Clearly this is not 
the case, (otherwise there would be no metals in na- 
ture), because quantum mechanically phases of the 
wavefunction add up coherently and the construc- 
tive interference of phases after the multiple scatter- 
ing over many sites results in a large final transi- 
tion amplitude. Thus the propagation through the 
periodic lattice. Introducing disorder on the periodic 
potential of the lattice destroys the constructive in- 
terference among the phases, which results in the 
phenomenon of localization. For this reason, Ander- 
son localization is a display of quantum mechanical 
multiscattering effects [42j. The averaged localization 
length of the system is obtained from the exponen- 
tial decay of the retarded Green's function and given 
by the ensemble average of the norm of the retarded 
Green's function G^ 1 , (7 =< jj >), by I =< lj >, 
1 = i < lMG-^xAWN >^ (i) ~ (2f) where 
localization lengths lj are related to jj by lj ~ — 1= 

and L is the size of the landscape sector, L ~ Nl p . 

Replacing the solution EqnJSJand Eqn|Sl (where for 
simplicity Sj is taken zero), back to WDW Eqn. (JBJ , 



we get Fj(a) ~ e ± V / ^' M .The wavefunction of the uni- 
verse solution is localized in the moduli 'superlattice' 
around some vacua Xj but it has an oscillatory be- 
haviour with respect to a. 

a) ~ - v L^ e ±l ^ a -^ 1 (10) 

As in p| we obtain the solution for Fj (a) in Eqn. © 
by using the Vilenkin boundary condition of choosing 
the outgoing modes at future infinity [2£|. The time 
parameter is obtained from Eqn. [? ] 

^/i~a = Hjt (11) 

where Hj is the expansion rate experienced by the 
local observers bound to the universe \tj. The os- 
cillatory behaviour needed for interpreting a as the 
time parameter is exhibited only for positive energies 
ij. The choice of outgoing modes boundary condition 
thus picks only wavefunctions localized around sites 
with positive energy ej > for 'giving birth' to a uni- 
verse. (AdS-)solutions with ij in the interval [-W, 0] 
are decaying solutions and thus physically irrelevant 
for this choice of boundary conditions as they can not 
give rise to a universe. That is, gravity (L g ) in Eqn.© 
breaks the 'particle-hole' symmetry e — > — e on the 'su- 
perlattice', namely:boundary conditions and L g dis- 
tinguish states of positive energy by exclu ding their 
symmetric negative energy partners|lOt 

112113. We 

briefly discuss at the end why this important symme- 
try for the behavior of the density of states near zero 
energy levels is closely related to the role of the bound- 
ary conditions on the selection of the most probable 
wavefunction of the universe and its A. 

Density of States and the Most Probable Uni- 
verse on non-SU SY Sector : We still have to an- 
swer the question: Which solution of Eqii.© is 
the most probable. At first this seems a hopeless 
task since energies are distributed randomly in the 
allowed energy band. However, the question can 
be addressed statistically by maximizing the den- 
sity of states (DoS) p{ti) ■ The single-particle av- 
eraged density of states can be obtained from the 
imaginary part of the advanced Green's function, 
ImG A {ej) ~ {e _^ 2+1 -i with poles at |e| = |Aj - i~fj\ 

through the expression p(e) = - < l\ImGji\N > 
or more explicitly in RMT from p(e) = ■h < 
TrS(e-H(x)) >„„= jL. f D(H x )P(H x )Im(G A ). For 
our simple 1-dimensional 'superlattice' this equation 
yields 



m * ^ (12) 

The dimensionality of the 'superlattice' does not 
affect the main inverse power law behavior of p on 
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energy but it does change the power of e, I in the de- 
nominator of Ean ll2l From Ean ll2l the maximum is 
around £j = e w 0, or Xj sa 7,-, Fig.l. Note that for 
a highly random system, the disorder strength can be 
as large as 7 ~ O(Vw). 

The most probable solution for the wavefunction of 
the universe V^o can now be found from the maximum 
of the density of states since p(e) provides the distri- 
bution probability of states. From Ean ll2l it can be 
seen that p(e) is peaked around ej « which leads 
to the conclusion that, with our approach, the most 
probable universe for the non-SU SY sector, is the one 
with a 'physical' cosmological constant J^-eff = e — 0. 
Although in Fig.l for the sake of generality, the plot 
of the DoS p(e) is shown for all energies, the rele- 
vant levels here consistent with the choice of Vilenkin 
boundary conditions are only the positive energy lev- 
els, e > 0. 

Let us define the terminology physical vs. bare 
here and distinguish between the 'physical' cosmo- 
logical constant e} = A e ff,j and the bare vacuum 
energy Xj — Ab are j. Our argument below is that 
the wavefunction of the universe experiences the self- 
gravitating properties of only A e / / = e since this is the 
relevant energy that enters in Ean llll for the expan- 
sion of that universe. Let us think of two 'observers': a 
'local observer' bound to the wavefunction of the uni- 
verse ^j[a, Cj], (while oblivious about the landscape 
vacua energies [Aj)], measures ej as his/her physical 
cosmological constant in his/her universe because it is 
only ej that determines the expansion rate Hj in Fj (a) 
from \JijOL = Ht\ but a ' superobserver' bound to 
the landscape superlattice will however observe all the 
vacua of the landscape and 'notice' that the wavefunc- 
tion $ j [a, €j ~ 0] is localized around some landscape 
vacuum with some large vacuum energy Xj , (thus the 
name 'bare vacuum energy' for Xj, although its per- 
turbed energy ej » 0). 

DoS is the most useful quantity for extracting statis- 
tical predictions about the landscape in this approach. 
Let us analyze briefly its main characteristics: The 
density of state is a maximum and does not diverge at 
e = due to the breaking of the 'particle-hole' sym- 
metry by gravity. If this symmetry were preserved 
the divergence of p(e) around e ~ would signal ei- 
ther delocalization of "0(0), I — 0(L = Nl p ) or the 
need to go to higher order corrections such as density- 



density correlations < 5p 2 >; We suspect that the 
density of states is a maximum around e = because 
of quantum effects like tuneling whereby the most lo- 
calized wavefunctions are the states with the lowest 
energy. Probably when the wavefunction gets local- 
ized on a site Xj with some large energy Xj, because 
of the short range interaction or tuneling, there is a 
large probability that this wavepacket will still sample 
the nearby vacua for awhile by tunneling towards the 
ones with lower and lower energies until it settles to 
the lowest 'physical' vacua energy level |44|:DoS falls 
off as power-law rather than exponential towards the 
tail end which would give a nonnegligible probability 
weight to the higher energy states. 

Discussion: String cosmology is still a field in 
making. A lot of effort has been devoted to it based 
on the expectation that the field will result in predic- 
tions and provide answers for some deep challenges 
facing cosmology. Progress has been made in address- 
ing puzzles like dark energyfTsL Il9|. pioneering work 
on string inflation [T2I Hp or even searching for ob- 
servational signatures 20] . However lately it was re- 
alized that progress achieved thus far soon may be- 
come a stumbling block:while rich in phenomenol- 
ogy, the field can provide not only one vacua solution 
for a universe like ours but many billions of billions 
of them, quite likely a whole landscape with 10 100 
vacua or morels H |H |H IH HI- The danger 
of losing predictability and thus becoming a nonfal- 
sifable theoryEH places string cosmology temporar- 
ily in an uneasy position from where even appealing 
to anthropic selection rules may seem an attractive 
possibility El El E3- 

The purpose of our approach, first introduced in 
, was to propose a new selection criteria for the 
landscape vacua, based on the dynamics of the wave- 
function of the universe propagating on the landscape 
background. We proposed to place quantum cosmol- 
ogy on the landscape background, thereby calculating 
the most probable wavefunction of the universe from 
solutions to the WDW eauationpl l32| .Although min- 
isuperspace is a limited aproach and issues like the 
role of boundary conditions and normal-ordering of 
operators are still under debate[2l IH |2l |^, l33l ? 
],this approach presently provides a dynamical quan- 
tum selection criterion for predicting the most proba- 
ble universe with vacuum energy A. Initially this selec- 
tion rule was applied to the SUSY sector of the land- 
scape and the most probable universes were 'stand- 
ing wave' solutions extended over the whole landscape 
SUSY sector, peaked around energies A ~ -j^, with 
N refering to the number of vacua[l| . 

The same analysis is applied in this work to 
the investigation of the minisuperspace restricted to 
the non-SUSY sector of the landscape and flat 3- 
geometries.In the absence of knowledge about the 
detailed structure of the non-SU SY sector, vacua 
energies were considered to be a stochastic vari- 
able, namely randomly drawn from the interval 
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[0, ±VF]. Wavefunction solutions found from WDW 
equation, Ean ll2l exhibit the well-known phenomenon 
of Anderson localization characteristic of disordered 
systems. Localization of the wavepacket ensures that 
coherence of the wavefunction of the universe is main- 
tained over large time-scales. Due to the stochastic 
distribution of landscape vacua energies Xi, i = 1, ..N 
we can not predict on which vacua of the superlat- 
tice the wavefunction \&(x, a) will be localized or the 
exact vacuum energy it will find there. However, 
we can make statistical predictions for the distribu- 
tion probability of the physical cosmological constants 
A e //,j = ij and the most probable wavefunction of 
the universe \tj[a, e?]. This information is extracted 
from the density of states p(ej). In the non-SUSY 
sector of the landscape our findings indicate that the 
most probable wavefunction of the universe solution 
selects states of zero energy and consequently of zero 
'physical' cosmological constant, (although the super- 
latice 'bare' vacua energies Xi where the most probable 
wavefunction localizes, may be as large as the disor- 
der strength W, Xj s» 0(7))- Statistically, the most 
probable 'physical' cosmological constant 'as seen by 
the ^Pj-bound local observer' through measurement of 
his /her expansion rate A e / / j = Hj = Wej, is the per- 
turbed energy |e 3 -| ~ |Aj±\/r| peaked around |ej| ~ 0. 
This behavior of DoS may result from the following: 
as the wavefunction gets localized at some vacua Xj 
with vacua energy Xj it can still sample nearby vacua 
by quantum mechanical tunneling to neighbor sites Xk 
that have lower vacuum energies A^. This interaction 
contained in Vi results in energy corrections from dis- 
order 7j. After 'moving around' to sample the neigh- 
bor vacua the localized wavepacket can finally settle 
to the site with the lowest energy. The probability of 
the reverse path is much smaller and thus does not 
dominate the DoS p(e). We suspect it is for this reason 
that the low energy states dominate the expression for 
p(e) . However it should be noted that the width of dis- 
tribution in the plot for the DoS, in Fig.l, is of order 
the disorder T. Therefore, depending on the strength 
of disorder, the probability of solutions for universes 
with nonzero cosmological constants may still be sig- 
nificant. The breaking of the 'particle-hole' e — ■> — e 
symmetry by gravity ensure that DoS does not diverge 
at e — > which guarantees the existence of localized 
solutions at zero energies [Tol ITU]. It is interesting to 
try to understand the relation between the existence 
of a localized solution around zero energies with time- 



reversability symmetry determined by the boundary 
conditions. The Vilenkin boundary condition we em- 
ployed here, clearly has a preferred direction of only 
outgoing modes at time infinity. A different choice 
of boundary conditions that would mix incoming and 
outgoing modes thereby preserving time-reversability, 
would likely result in a divergent DoS near zero en- 
ergies,which can be interpreted as the lack of solu- 
tions with zero physical cosmological constant from 
the probability view of DoS. The relation between DoS 
and time-reversability symmetry deserves careful ex- 
amination and will be reported elsewhere since it is 
beyond the scope of this paper. As an aside, from the 
expression for DoS, Eqn^l it should be noted that 
states with > €j < 7 still have a non-negligible prob- 
ability since the occupation of states at the tail end of 
DoS is only power-law suppressed. 

DoS plays the role of an order-parameter for the 
disordered 'superlattice'. Its peculiar behaviour near 
zero and critical energies can signal phase transitions 
and critical phenomena. Therefore it is a very impor- 
tant quantity for exploring the rich phenomenology of 
the landscape. Many issues remain to be addressed 
yet. The limited degrees of freedom for the minisu- 
perspace we considered in this work do not allow ad- 
dressing problems like the coupling of the standard 
model to landscape moduli and its backreaction effect 
on the expansion rate Hi, the DoS and the most prob- 
able wavefunction, or the later evolution of a universe 
born out of the wavefunction localized around some 
landscape vacua. Future directions for the applica- 
tion of our proposal have to involve an extension of 
the minisuperspace to more degrees of freedom needed 
for addressing more realistic scenarios, for example an 
extension to non-flat 3 geometries or to universes with 
the standard model content in them. The latter, cur- 
rently under study, can be used to relate Higgs to A 
hierarchies through universal scaling and critical phe- 
nomena. We report our results about the latter in a 
separate publication. Despite the current limitations 
of quantum cosmology, we believe it is important and 
a step forward to provide an alternative proposal to 
the anthropic principle that dynamically selects the 
landscape vacuum in which our universe resides, as 
shown in this work. 
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